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FOREWORD 


In Part I of this report (Schwiderski, 1978a), a combined hydrodynamical-empirical method 
was introduced to compute numerically harmonic partial tides in the world oceans with an ac- 
curacy of better than 5 cm, which is needed in various military and civil applications of today. In 
this report, the computed semidiurnal principal solar tide (S, ) is displayed in an atlas of tabulated 
tidal charts and plotted corange and cotidal maps. 
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ABSTRACT 


“In Part 1 ¢(Schwiderski, 19784) of this report, a unique hydrodynamical interpolation tech- 
nique was introduced, extensively tested, and evaluated in order to compute partial global ocean 
tides in great detail and with a high degree of accuracy. This novel method has been applied to 
construct the semidiurnal principal solar (S,) ocean tide with a relative accuracy of better than 
Semiunywhere in the open oceans. The resulting tidal amplitudes and phases are tabulated on a 
t° x 1° grid system in an atlas of 42° x 71° overlapping charts covering the whole oceanic globe. A 
corresponding atlas of global corange and cotidal maps is included to provide the reader with a 
quick general overview of the major tidal phenomena. The specifying hydrodynamical parameters 
of the model are listed along with quoted sources of empirical tide data, and significant tidal 
features are explained and discussed. The S, ocean tide is found to resemble closely the 
corresponding lunur M, tide presented in Part I] of this report.- | | : 
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1. INTRODUCTION 


Part [| of this report (Schwiderski, 1978a) introduced a unique combination of hydrody- 
namical and empirical methods to model detailed ocean tides with a relative component accuracy 
of better than 5 cm anywhere in the open oceans. This enormous accuracy is well above minimum 
requirements set by. for instance, the National Aeronautics and Space Administration (NASA) 
and the Department of Defense (DoD) — to map the geoid at sea by satellite altimetry to within 
i0 cm. The following features of this unique hydrodynamical interpolation model made the 
achievement of this accuracy possible. 


a. A spherically graded 1° x 1° grid system is set up in connection with a corresponding 
1° x 1° bathymetry to assure a sufficient resolution of all important tidal phenomena. 


b. The bathymetry of the gridwise. simply connected ocean basin is hydrodynamically 
defined (Schwiderski, 1978c) by appropriate modifications of earlier realistic depth data col- 
lections. The hydrodynamical redefinition was needed in order to model the well-known strong 
distortion and retardation effects of shallow continental shelves, narrow ocean ridges or island 
chains, and other significant bottom irregularities. 


c. The Boussinesq substitution of the turbulent Reynolds stresses is applied in the form of 
eddy dissipation with a novel physically meaningful eddy viscosity that depends linearly on the 
lateral grid-cell area and, hence, directly on the ocean depth. 


d. The linear law of bottom friction is introduced with a bottom-friction coefficient depend- 
ing linearly on the bottom grid-cell area which is independent of the ocean depth. In boundary 
cells, the otherwise constant friction coefficient is subjected to an indirect cellwise adjustment in 
order to permit a consistent hydrodynamical interpolation (see h., below) of empirical tide data 
known from tide gauge stations at continental shores, islands, or other shallow-ocean bottom 
irregularities. 


e. The effects of the terrestrial tide and the oceanic tidal load are included as simple second- 
order approximations in the sense of Love and Accad and Pekeris (1978). 


f. The Hansen-Zahel (Zahel. 1970 and 1977; Estes, 1977) finite differencing technique ts 
modified by a new differencing scheme in time which improved decay, dispersion, and stability 
characteristics of the numerical procedure and facilitates the simple indirect adjustment of the 
bottom-friction coefficient in the hydrodynamical interpolation technique (see d. and h.). 


vy. At land-ocean cell walls. the conditions of no-flow across and free-slip along the bounda- 
ries ure enforced. The no-flow condition is subsequently relaxed by allowing controlled periodic 
inflows and outflows over the mathematically assumed boundaries. This allowance redetines 
indirectly more realistic shorelines in order to further improve the consistency of the hydro- 
dynamical interpolation of empirical data (see d. and h.). 
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h. A unique hydrodynamical interpolation technique is introduced which incorporates into 
the theoretical model empirical tidal constants collected from over 2 000 tide-gauge stations 
around the world in a hydrodynamically consistent fashion (see d., f., and g.. above). 


'. A new higher-order approximation of arctic ocean tides is used, that is described in 
Schwiderski (1981). 


With these features, the new model was successfully applied to chart the semidiurnal principal 
lunar (M3) ocean tide with the desired accuracy. The technique and accuracy of the mode] were 
extensively described and discussed in Part I of this report as well as in subsequent publications 
and symposia presentations by the author (Schwiderski 1978a, b; 1979a. b, c. d, e; and 1980). 


The same hydrodynamical interpolation technique has been applied to chart the semidiurnal 
principal solar (S,) ocean tide with the same relative accuracy as M,. Again, it must be emphasized 
that the enormous accuracy achieved over all open ocean regions diminishes somewhat near 
coastal areas where Known empirical data are marginal in quantity and/or quality. 


A complete listing of all sources of empirical ocean tide data, which were interpolated into 
the S, tidal charts, is presented in Appendix A. In the meantime, Section 2 of this report lists 
the significant hydrodynamical input parameters that specified the constructed S, ocean tide. 


The major features of the global S, tide are discussed in Section 3. A complete numerical display 


is presented in Appendix A where all tidal amplitudes and phases are gridwise tabulated in map- 
like charts. Corange (equi-amplitude) and cotidal (equi-phase) maps of the S, ocean tide are 
plotted in Appendix B. 
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2. S, OCEAN-TIDE PARAMETERS 


The astronomical semidiurnal principal solar (S,) equilibrium tide 7 (or tide-generating 
potential Gy: see Schwiderski, 1978a) at the geographical point (A. @) and instant (7) is determined 
Hy 


n = A cos? @ cos(or +X + 2A) (1) 
where 
G = 9.8) m/sec? earth gravity acceleration 


longitude (east in rad) 
latitude (north in rad) 


— 
— 


A 
7) 
f = universal standard time of day (in sec) 

K = 0.412 3841 m=S, equilibrium tide amplitude 
o = 1.45444 10 4 sec | =S, tide frequency 

x 


= 0=S, ustronomical argument (in rad) 
The corresponding instantaneous ocean partial tide (Schwiderski, 19784) is determined by 
f= Ecostart+X- 4). ‘ae 
where the focal harmonte constants 
£ = &A.6)=S, ocean tide amplitude (in m) 
and 
5 = 5(A. 6) = 8, ocean tide Greenwigh phase (in rad) 
must be determined, say. by linear interpolation in the tidal charts of Appendix A. 


A simple second-order approximation in the sense of Love and Accad and Pekeris (see Part. 
Schwiderski. 1O78a. 1979¢, and 1980: and Accad and Pekeris, 1978) yields 


f= 0.0127 and 6° = - 0.0607, 


ne. the corresponding terrestrial tide (© and the earth dip C°? (yielding) under the oceanic tidal 
load €. respectively. A more elaborate and probably slightly more accurate earth dip oO May he 
computed by using Farrell’s Green function (see Farrell, (972 and 1973. and Schwiderski. 1980). 
In linear superposition. one finds the corresponding instantancous geocentric partial Sy tide: 
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A detailed deseription of the hydrodynamical-empirical medel to compute the ocean tidal 
amplitudes & and phases 6 (listed in Appendix A) was given in Schwiderski (19784, 1979¢, d, and 
1980), In particular. afl model input parameters such as the dimensionless eddy coefficient € (q's. 
[O03 and [23). the bottom-friction parameter > CEq’s. 4a and b), and the differencing parameters 
K and kK (Eq’s. 64 and 72) were all specified in Schwiderski (1978a) (referenced equations). These 
parameters were determined for M, by extensive trial-and-error computations and remained 


unchanged tor the construction of S5. 


In the computation of the S, tide model, the following mode-dependent parameters were 
used (see referenced equations in Schwiderski, 19784): 
a. The time step Ar (Eq’s. 64, 123): 


Ar = 180.000 see (S) 


b. The hydrodynaiical interpolation control limits, A,.A,, and A; (Eq’s. 88, 89, 94. 97, 
and 99) rematned the sane as for M,:; 


hk, = 0.045, 4, =0.045. 4, =0.5. (6) 
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( 3. S; OCLAN-TIDE FEATURES 
Phe entire constructed Ss ocean tide is gridwise displayed in map-like amplitude and phase 
tables mn Appendix AL Phe 42° \ 71° charts cover the whole globe north of colutitude 169° 
(Antartica) in three zones: a northern zone N from O° to 71° colatitude. a middle zone M from 
48° to TES colutitude. and a southern zone S from 98” to 168° colautitude. The overlapping 
voourapliical areas of the Gdal charts have been chosen to provide a worldwide coverage for special | 
dppheations and to allow the reader to scan the large amplitude and phase charts together in order 
to evatite then quality and visualize the important tidal features. In addition. a generally super: ) 
Hho overview ob sore tidal features can be recognized by inspecting the more schematically | 
plotted corange and cotidabiaps provided in Appendix B. 


boran casy evaluation of the tidal charts in Appendix A, all hydrodynamically interpolated 
cmpirtcal tidal amplitudes and phases have been visibly murked by subbars tor all shore data and 
subbrackets: tor all nearshore deep-sea input constants. Furthermore. the charts display the 
approsimnate locations of distant off-shore deep-sea stations by subtildes under the computed 
amplitude and phase data, Phe corresponding empirical data, which were excluded from hydro- 
dynamical interpolation (see Seet. | and Schwiderski. 19784. 1979d. and 1980), are listed and 
compared with the modeled data Fables bo 2. and 3. Finally. the approximate geographical 
locations of the important amphidromic points of zero amplitudes are marked by a circled 


The tidal charts and maps permit the viewer to follow the tidal waves. that is the high water 
fronts: Cerests), in forward Cor backward) direction. for tastance. on their rotation around the 
amphidromie points. In the tidal phase charts of Appendix A. it is best to start from the pro- 
minently visible 07 = 360° or 100° cotidal lines. Since the Greenwich phases specify the time lags 
(in degrees: 300 = To hour) of the tidal crests relative to the cresting time of the corresponding 
equilibrium tide along Greenwich meridian, one gathers a vivid impression of the significant global 
and local tidal phenomena, 


By following the tidalwaves on their periodic rotations, One finds these waves passing through 
the specmlly marked stations in empirically correct time and with the correct height. In tact. all 
over the globe over 2 000 tidal phases and 2 000 amplitudes are coherently integrated. This is 
particularly tapressive for the charts of the Pacific Ocean. where the empirical data from so many 
clustered and scattered island stations fit smoothly into the surrounding computed tides. From the 
smoothness features of erratically interpolated tidal data (see Parts Pand HH). one concludes that 
this results notan artifact of the interpolation applied but constitutes a vivid manifestation of the 


excellent compatibility of both the empirical and hydrodynamical procedures combined. 


On the basis of Chis observation, it can again (see Schwiderski, 197&a. be 1979a, bo dee: 
and TOSO) be estimated that the Sy tidal charts permit a tide prediction with a uniform accuracy 
relate to VE, of better than 5 car anywhere in the open oceans. Naturally. near rough ocean basin 
rohets Geog. Aretig and: Antarctic shores), where empirical tide (and depth) data are margmal tn 


quality and quontiey. a somewhat lesser dccuraey must be expected. The estimated aecuracy ot 
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the computed S, tide ts, of course, fully validated by all 32 empirical tide data from distant off- 
| shore deep-sea Ude gauge stations. which are listed along with the computed duta in Tables |, 2, 
und 3. The differences (not necessarily errors) range from 0 to 3 cm in amplitudes and 0° to 6° 
(12 minutes) in phases (with one 13° exception) and thus verify the estimated prediction accuracy. 


From the tidal charts and maps in Appendixes A and B, one concludes the following major 
features of the S, ocean tide, which are very similar to M2 (see Schwiderski 1979e): 


Table 1. North Atlantic Ocean Deep-Sea Empirical and Modeled S, Tides 


LONGW LATN:EMP£& MODE At. EMP5 MODS Ad IAPSONR’ SOURCES 

13°51’ 58°16’ | 32 33 +{ 204 207 $3 1.1.37 C 

24°43' 62°50’ 5 35 35 Ges 39 22043 1.1.29 C 

28°46’ §=60°12’ | 27 2°80 +1 $222 225 +3 1.1.30 C 

2938 SFO 19 i ene 222, 1.1.31 C 

30°10’ = 53°39" 3 12-1) 198 201. -+3 eee a 

25°06’) = 53°31’ 18 18 0. 183 S70 +4 1.1.33 C 

20°00’ ~=—- 3°39" 25 +] 7? 179 +2 1.1.34 C 

28°11’ 48°45" | 10 1] +] 14] l42 +] 1.1.38 C 

78°09’ = 45°? ]’ | J? 13 +] 112 114 +2 1.1.39 C 

27°57’ = 44°25’ | 15 17 +) g? 94 +2 1.1.40 C 

20°05’ =. 37°09’ 233 24 +1 83 88 +5 11.41 C 

14°15' =36°41' «30 30 Qi 9? 89 3 ie es. C 

5°38 32°47 9 8 +] 5 27 % -'iI i 3 C.M 

1635. 20°26") 8 8 Q'° 23 25 42 Lor CLP 

76°48’ 28°27" 7 7 Q- 32 34.0—Ca+D? 1.2.15 C 

76°47' 28°01’ 7 7 O: 4 35 OO C 

67°32" 28°14 7 7 0, 2322 3300] ee aes | 
69°45' = 28°08" - 7 7 QO, 3 33 +> ae CZ : 
69°40’ 27°59" 7 6 +f. 30 34.00 +4 i a CZ | 
69°40’ 27°58") 7 6 —-| i; 33 34+} 1.2. 7 C.Z ! 
69°20’ 26°28" | 6 Q! 34 350 #2) 4.2.10 CZ | 
69°19" 96°27’ 7 6 —| 3} 35 +4 E25 C.2 


= Aniplitudes (cm) 
= Greenwich Phases (dey) 
APSO = Int. Assoc. for the Phys. Sci. of the Oveans , 
€ = Cartwright et al. (1979) 
| M = Mofjeld (1975) 
| P = Pearson (1975) 
7 = Zetleret al. (1975) 
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Table 2. Northeastern Pacific Ocean Deep-Sea Empirical and Modeled S, Tides 


LOGW LATN , EMPE MODE A¢ | EMPS MODS Ad I[APSONR SOURCES 
ee a, ee ee —4>- Be Sects hd a He 

144922’ 56°08" | 32 33 “e 315 ce ue Mal ¢ 

135°38' = 53°19’ | 30 33 +3: 30) 303. +2 21.16 C 

132°47' 49°35' 29 28 -} +) 283 2870 +4 21.15 ( 

145°00' 34°00' | 12 — 323 

145°00' 34°00':—- {2 a4 323 

124296" 274s". 12-2" 116 120 +4 13 CoM 

129°O1' 24°47" |? 9 -3! 95 oy, 22 1.40 CoM 


g = aeipitudes (cm) 

6 = Greenwich Phases (deg) 

[APSO = Int. Assoc. for the Phys. Sci. of the Oceans 
C = Cartwirght et al. (1979) 

M = Munk et. al. (1970) 


Table 3. Southeast Indian Ocean Deep-Sea Empirical and Modeled S, Tides 


LONG E LATS | EMPE MODE st ; EMP5 MODS AS IAPSONR SOURCES 


= a eA seh SE a re et las “a $2 5. ote we 
132°01' 37°01’) 16 17 +1: 100 104s +4 a CUS 
132°09' =©50°02’ | 9 9 QO: 108 114 +6 2 ae: CIS 
132°07' 60°01’ ; 10 9 ~1 | 120 133 #13 7 Clas 


§= Amplitudes (cm) 

& = Greenwich Phases (deg) 

JAPSO = Int. Assoc. for the Phys. Sci. of the Oceans 
é = Cartwright et al. (1979) 

[S = Jrish and Snodgrass (1972) 


a. Strong, Weak, and Scie Amphidromes, These exist amphidromie rotations around 
so-called “strong, “Weak,” “incomplete” amphidromes. The fully developed strong amphi- 
dromes display almost cee INcreasing amplitudes and phases im radial or tangential directions. 
respectively. An example of the strong amphidrome type may be seen mn the North Atlantic 
(Table 9N). These amphidromes are distinctly isolated and independent of adjacent amphidronn, 
rotations. 


Weak amphidromes are characterized by distinctly unevenly growing amplitudes and phases. 
These amphidromes occur mostly in pairs with a channel of low amplitudes and rapidly varying 
phases in between. A typical example can be seen in the central bast Pacific (rable 7M). 
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binally, there exist areas with distinetly low amplitude tides and rapidly varying phases that 
seem to converge toward one €or more) points (incomplete amphidromes). k-xamples of this type 
may be seen in the South Atlantic (Fable 9S) and in the center of the Caroline Islands (Table 4M ). 


b. Compatible and Incompatible Amphidrome Pairs, There exist pairs of amphidromes with 
“contpattbly” (oppositely ) or “incompatibly” (parallely ) rotating tidal waves. An example of the 
first Kind, located in the North Pacific. has a weak clockwise amphidrome in the center and a 
counterclockwise amphidrome in the east. An example of the second kind exists in the North 
Atlantic. This pair has a strong counterclockwise amphidrome in the center and a neighboring one 
In the eastern Caribbean Sea. 


As a result of incompatible amphidromes. one finds lines in cotidal charts along which 
progressing tidal waves become standing and reverse their directions. On those wave reversal 
curves, coudal lines are cuspate. Another interesting consequence of incompatible amphidromes Is 
the progressing tidal wave that follows along the whole eastern coast of the South and North 
Atlantic from South Africa in the south to Norway in the north. This tidal wave is obviously 
driven by the combined amphidromic systems existing in the South Atlantic. in the eastern 
Caribbean Sea, in the North Atlantic, and north of the Faeroe Islands, all of which rotate counter- 
clockwise. 


C. Relief Distortions and Retardations of Tides. As can be seen in all amplitude and phase 
charts of Appendix A, tidal waves crossing shallow continental shelves (see. e.g.. the western 
Atlantic coast), narrow ocean ridges. or island chains (see. e.g.. Caribbean, Hawaiian. Aleutian, or 
Marianas) suffer strong distortions and retardations. In fact, phase retardations of more than 30° 


and even 60° occur across ocean ridges. Moreover in certain coastal corners (e.g., in the North 


Sea. Table IN: in the Bering Sea, Table SN; on the Patagonian Shelf, Table 8S: and northwest of 


Australia. Table 4M). one finds oceanic tidal waves shedding small amphidromic rotations. As was 
pointed out in Schwiderski (1978a, 1979d, e. and 1980), those rough tidal phenomena have long 
been known (see e.g., Harris. 1904; Defant, 1961: and Luther and Wunsh. 1975). The satisfactory 
resolution of those distortions played a major role in the successful charting of ocean tides. 
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4. CONCLUSIONS 


The hydrodynamical interpolation technique hus been applied to construct the semidiurnal 
principal solur tide (S,) with a relative accuracy of better than 5 cm anywhere in the open oceans. 
Phe constructed tide ts displayed by tabulated charts in Appendix A and by corange and cotidal 
Maps in Appendix B. The major features of the S, tide ure discussed in Section 3. A comparison 
with the earlier computed M, model reveals close similarities. 
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APPENDIX A 


ATLAS OF 1° x 1° S, OCEAN TIDE AMPLITUDE 
AND PiAASE Cd ARTS FOR 42° x 71° AREAS 


f. GUIDE TO TIDAL CHARTS 


M = m: Longitude Number 
N = n: Colatitude Number 
Nn = (m ~ 0.5)°: Geographical Longitude East 
0, = (n- 0.5)°: Geographical Colatitude 
Emin = EA. 8): Amplitude (in cm) 
bmn = S(Ag, . 8): Greenwich Phases (in deg.; 30° = 1 hr) 
@) = Amphidromic Points 
= Subbars Mark Empirical Input Data at Shore Stations 
4 = Subbrackets Mark Empirical Input Data at Near-Shore Deep-Sea Stations 
a = Subtildes Mark Approximately Distant Offshore Deep-Sea Stations with Excluded 


Empirical Tide Data Listed in Tables 1. 2. and 3 
-. SOURCES OF EMPIRICAL TIDE DATA 
Publications: 


National Ocean Survey (1942). British Admiralty (1977). International dydrographic Bureau 
(1978). Defant (l9GE). Miyazaki et al. (1967), Nowroozi et al. (1969). Munk et al. (1970). Zahet 
(19°70), Irsh et al. (1971). Trish and Snodgrass (1972), Nowroozi (1972). Luther and Wunsh 
(197S). Motyeld (1975). Pearson (1975), Zetler et al. (1975). Cartwright et al. (1979), and Pugh 
(L979). 


Private Communications: 


ID. ©. Simpson (1977). National Ocean Survey. Rockville, Maryland: S. K. Gill and D. L. Porter 
(1978). National Ocean Survey. Rockville. Maryland; K. Wyrtki (1978). University of Hawaii. 
Honolulu. dawan. and D. b. Cartwright and D. Pugh (1978). Institute of Oceanographic Sciences. 


Bidston Observatory. UK 
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ATLAS OF CORANGE AND COTIDAL MAPS 
. OF THE S, OCEAN TIDE 


Amplitudes & of corange lines in cm. 


Greenwich phases 6 of cotidal tines in 15, 30, 45, 60. 75. 90. 105, 120, 135. 150. 165, 
180, 195. 210, 225, 240, 255. 270, 285, 300. 315, 330. 345. 360 = 0° where 30° = 1 hour. 
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